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Q\ . Abstract 

CO ' We demonstrate that the generalized Gell-Mann-Low theorem permits for a sys- 

tematic expansion around the nonrelativistic limit when applied to bound states in the 
, Wick-Cutkosky model, Yukawa theory, and QED (in Coulomb gauge). We apply this 

expansion to obtain new results for the fine and hyperfine structure of bound states 
in the cases of the Wick-Cutkosky model and Yukawa theory, and reproduce correctly 
Ph, the fine and hyperfine structure of hydrogenic systems. 

CD . 

1 Introduction 

> ■ 

^ ■ In this paper we present a method to extract the fine and hyperfine structure, presently to 
lowest order, of bound systems appearing in general relativistic quantum field theories in 
a systematic, consistent and transparent way. It is based on the application of a general- 
ization of the Gell-Mann-Low theorem pQ which generates an effective Hamiltonian called 
the Bloch- Wilson Hamiltonian, for the dynamics of the two constituents. The Bloch- Wilson 
Hamiltonian consists of the relativistic kinetic energy of the constituents and an effective po- 
tential that can be obtained explicitly in a perturbative expansion in powers of the coupling 
constant, derived from the Dyson series for the adiabatic evolution operator. At every per- 
turbative order, the effective potential can be expanded in powers of the relative momentum, 
which, together with the corresponding expansion of the relativistic kinetic energy, leads to 
a systematic expansion around the nonrelativistic limit. 

We will demonstrate this procedure in three different theories, a scalar model with cu- 
bic couplings (the "Wick-Cutkosky model"), Yukawa theory, and quantum electrodynamics 
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(QED). The Bloch- Wilson Hamiltonian for two of these theories has been carefully derived 
to lowest nontrivial order (in the expansion in powers of the coupling constant) and used for 
numerical bound state calculations before (21 El • The Bloch- Wilson Hamiltonian for QED 
in the Coulomb gauge has been discussed very briefly in Ref. [3], and will be treated in 
much more detail in a future publication. In each case, we will take the exchanged boson 
to be massless in order to have analytic expressions in the nonrelativistic limit and for the 
lowest-order fine and hyperfine structure. 

This procedure should be compared with the historical calculation of the complete fine 
and hyperfine structure in hydrogenic systems. Ref. [5 j gives a good overview of these cal- 
culations which are based on the Breit equation [S| and have not been completed until 1951. 
A more systematic treatment is provided by the Bethe-Salpeter equation [7J |K| in the per- 
turbative expansion devised by Salpeter [Hj. However, Salpeter's expansion begins at zero 
order with an instantaneous interaction of the constituents which is natural in QED, par- 
ticularly in the Coulomb gauge, but not in the case of spinless boson exchange as in the 
Wick-Cutkosky model and Yukawa theory. In fact, the analogue in the Bethe-Salpeter ap- 
proach to the procedure outlined above would be to consider the lowest-order approximation 
to the Bethe-Salpeter kernel, the so-called ladder approximation, and expand the equation 
in this approximation around the nonrelativistic limit. However, even in the simplest case 
of the Wick-Cutkosky model [TU] this procedure leads to a "curious" term of the unexpected 
order a 3 In a jTTj. The full fine structure in this model and in Yukawa theory have, to our 
knowledge, never been calculated. 

In the next two sections of this paper, we will present the Bloch- Wilson Hamiltonians for 
the field theories mentioned above. In fact, we will discuss two different types of effective 
Hamiltonians, in Section 2 the Bloch- Wilson Hamiltonian used before in Refs. (21 El H] which 
generalizes in a natural way Bloch's formulation of degenerate perturbation theory J2|, and 
in Section 3 a hermitian version used before by Wilson ^3] and, independently, by Gari et al. 
[T4] and by Kriiger and Glockle building on earlier work by Okubo 16 . As will become 
clear in Section 4, the two types of effective Hamiltonians differ by an antihermitian term in 
the expansion around the nonrelativistic limit which is present in the first type of effective 
Hamiltonians (except for the case of QED) but obviously not in the second. In Section 4, 
we will describe how to calculate the fine and hyperfine structure to lowest order starting 
from these effective Hamiltonians, relegating the more technical aspects of the calculations 
to Appendix El In particular, we will use second-order time- independent perturbation the- 
ory to calculate the contributions originating from the antihermitian term. An alternative 
method based on a Foldy-Wouthuysen-type transformation from a nonhermitian to a hermi- 
tian Hamiltonian is detailed in Appendix El We will also comment on the possible physical 
significance of the antihermitian term. In Section 5, we discuss the results for the fine and 
hyperfine structure in the different theories and compare in particular with the numerical 
solutions (21 E] for the Bloch- Wilson Hamiltonian (to lowest nontrivial order in the expansion 
in powers of the coupling constant) without the expansion around the nonrelativistic limit, 
and with the well-known fine and hyperfine structure in hydrogenic systems. 
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2 The Bloch- Wilson Hamiltonian 



We will begin by stating the generalization of the Gell-Mann-Low theorem proved in Ref. 
PP. To this end, some general notations have to be introduced first: the full Hamiltonian H 
of the field theory under consideration is decomposed into a "free" Hamiltonian H (typically 
the one describing free particles) and an interaction Hamiltonian Hi, H = Hq + H\. The 
adiabatic evolution operator U e maps a state \<f>) in the Fock space T in which H Q acts, to 
the state that results from evolving \<p) with the Schrodinger equation corresponding to 
the adiabatic Hamiltonian 

H t {t)=H, + e-^H ll e>0, (1) 
from £ — ► — oo to £ = 0. The adiabatic evolution operator has the perturbative expansion 

U e = Yl —T / d *i ■ ■ ■ / 4e- £(|tl,+ - +|t " l) T[F 1 (t 1 ) ■ ■■H l {t n )) , (2) 

the well-known Dyson series, where 

#x(£) = e iHot H 1 e- iHot . (3) 

For the generalized Gell-Mann-Low theorem, a linear ifo-i nv ariant subspace Qq C JF, 
H Q C f2 , is fixed and the corresponding orthogonal projection operator P : ^ — ► 
introduced. The theorem asserts that, if the operator 

U B = lim U^PoU.Po)- 1 : fi ^ ^b(^o) = , (4) 

e— >0 

the Bloch- Wilson operator, is well-defined in Qq, then the image subspace Q is invariant 
under H, HQ C fi. By its definition, Ub is "normalized" to 

Pot/B = lno = ^olo , (5) 

consequently the inverse of Ub is Po\q- This normalization naturally generalizes the one 
used in the original Gell-Mann-Low theorem jH] and corresponds to the one used by Bloch 
in his formulation of degenerate perturbation theory ^2]- Note that we have changed the 
notation for the Bloch- Wilson operator compared to our previous works to emphasize this 
fact, for reasons to become clear later. The existence of the operator Ub and the assertion 
of the theorem are to be understood order by order in the corresponding formal power series 
originating from Eq. (J2*|). as is the case for the original Gell-Mann-Low theorem. 

The fact that Q is if -invariant is equivalent to the diagonalizability of H in Q, so that part 
of the eigenvalue problem of H can be solved in Q. It is convenient to similarity transform 
the eigenvalue problem back from Q to the usually more manageable subspace Qq. Hence 
we define the effective or Bloch- Wilson Hamiltonian 

H B = P HU B : fio -> • (6) 
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We can use Eq. (J5J) and PqHq = H P as a consequence of the if -invariance of fl , to bring 
Hb into the form 

H B = P H U B + Po^Ub 

= H P + PqH.Ub . (7) 

The eigenvalue problem of Hb in f2 is equivalent to the one of H restricted to Q with the 
identical eigenvalues. The corresponding eigenstates are mapped into each other by U B and 
P . 

For the following bound state calculations we will take Q to be the subspace of all states 
of the constituents as free particles (under Ho), usually in a momentum eigenstate basis. In 
this case, the first part H P in Eq. (JJJ) just represents the kinetic energy of the constituents, 
consequently the second part defines an effective potential for their interaction (although it 
also contains radiative corrections to the vacuum energy and the masses of the constituents) 
given in terms of a power series in the coupling constant by inserting the Dyson series (J2J). 
The operator P projects out from the true if-eigenstates in O the component with the 
particle numbers corresponding to the constituents, hence we can think of the corresponding 
ife-eigenstates as wave functions of the constituents, the higher Fock space components 
(with higher particle numbers) being generated by the application of Ub- For bound states, 
it is these constituent wave functions that have to be normalizable, while the corresponding 
if -eigenstates will usually not be normalizable with respect to the standard Fock space scalar 
product. 

In the following, we will calculate Hb to lowest nontrivial order for two-particle subspaces. 
In all field theories to be considered, the interaction Hamiltonian H\ changes the particle 
number so that 

P H 1 (t)P = 0. (8) 

In this case, the lowest-order nontrivial contribution to Hb arises from the first-order term 
in the expansion (J2J), and we can write 

H B = H P - i f dte-^PoH^H^Po + 0{Ht) , (9) 

J — oo 

where the limit e — > is understood. 



2.1 Wick-Cutkosky model 

We will now present the Bloch- Wilson Hamiltonian Hb for three different theories. We start 
with the "Wick-Cutkosky model", a field theory with three scalar fields 4>a, 4>b, and <p, the 
latter taken to be massless, with interaction Hamiltonian 



Hi = g d 3 x : 



it 



x + 



V( x ) 



(10) 
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We consider states with one A and one B scalar. The Bloch- Wilson Hamiltonian for the 
corresponding subspace was calculated in Ref. j2], with matrix elements 



(p A , Pb\H b \p' a , p' B ) = (^m 2 A + p 2 A + y/wl + plj (2n) 3 5(p A - p a )(2tt) 3 5(pb - p^) 
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2E£ 2EZ 2E A , 2E B , 

PA Pb p' A p' B 



X 2\ PA - p> A \ { E£ A + \p A - p' A \ - E$ a + EB b + |p B - p ^| - E* b 

x (2nf5{p A + p B -p' A -p' B ) (11) 
in a momentum eigenstate basis with the nonrelativistic normalization 

(Pa,Pb\p' a ,p' b ) = (2tt) 3 5( P a - p A )(27r) 3 ,5(pB - p' B ) . (12) 
We have introduced the shorthands 



K' = \j m %B + V 2 (13) 

for the relativistic kinetic energies in Eq. (fTTj) . To arrive at the form the vacuum energy 
(including its lowest-order radiative corrections) was subtracted, and radiative corrections 
to the kinetic energy have been absorbed in a renormalization of the masses. The renormal- 
ization procedure was analyzed in detail in Ref. |3] from the present not manifestly covariant 
point of view, and the appearing one- and two-loop expressions have been shown to coincide 
with the ones of usual covariant perturbation theory. The renormalized masses are denoted 
here as m A>B - Making use of the overall momentum conservation Pa + Pb = p' A + p'bi the 
dynamics can be reduced to the center-of-mass system p^ + p# = 0. One then has for the 
corresponding Schrodinger equation, with p = p^ = — Pb, 

g 2 f d 3 p' 1 



m 2 A + p 2 + J m 2 B + p 2 J 0(pj 



2E£2EBJ ( 27r ) 3 J2E52E5 



'p -~p y -~ P ' -~ P ' 

\ I I 1 



X 



+ „».,.. 31 55 = ^(P) » ( 14 ) 



2|p-p'| ^+|p-p'|-^, ££ + |p-p'|-.E* 
with the wave function 0(p) defined as 

(Pa, Pb|0> = 0(pa)(27t) 3 5( P a + p B ) • (15) 
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2.2 Yukawa theory 

As our second example, we consider Yukawa theory, consisting of two Dirac fields ip A , i>B, a 
massless scalar field (p, and the interaction Hamiltonian 



Hi = g J d 3 x : [ip A (x)ip A (x) + V , b(x)^b(x)] cp(x) : 



(16) 



The Bloch- Wilson Hamiltonian for bound states of one A and one B fermion was calculated 
in Ref. j^j. The matrix elements are 

(Pa, r; Pb, s\H B \p' A , r'\ p' B , s') 



m A + P\ 



mi 



b + Pb W^Pa - p' A )S r A27r) 3 5(p B - p' B )8 s 



2E A 2EZ 2E A , 2E B , 

PA Pb p' A p' B 



2\p A - p' A \ \ + \p A - P^| - E* E* + \p B - p' B \ - E* 



x [u A (p A , r) u A (p' A , r')] [u B (p B , s) u B (p' B , s')} (2tv) 3 5(p a + p B - p' A - p' B ) . (17) 

The parameters r,s with possible values 1,2 in the momentum eigenstates \p A ,r;p B , s) 
describe the spin orientations of fermions A and B, respectively. The Dirac spinors u A (p,r) 
are related to Pauli spinors Xr via 



u A {P,r) 



with the Pauli spinors normalized to 



Ep + m A 



I Xr 

p ■ a 



\ 



Xr 



{11 



\ -E'p + m A ) 



X^rXs $rs ■ 



(19) 



Note that Eqs. (fTTjl and (fT7|) only differ by the Dirac spinor products (and the additional 
spin indices). 

The remarks after Eq. (fTTjl apply equally to Eq. (j!7jl . In the center-of-mass system, we 
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can write the corresponding Schrodinger equation conveniently in Pauli spinor form as 



™ A + P + \ ™b + P 0(P 



., ,E£ + m A EZ + m B r d 3 p' E$ + m A E$ + m B 
■ f ] 1 ~^2E£ 2Ef~ J (2^F V 2Ej~ 



+ 



2|p - p'| \E£ + |p - p'| - E$ EB + \p- p'| - E$ 

xfi- v A aA p '' aA )(i- p R aB V ' B CTB )m = Em, (20) 

^ E£ + m A E$ + m A J\ E$ + m B E* + m B J nP J ' 1 ; 

where the spinorial wave function 0(p) in the center-of-mass system is defined as 

<t>(p A )(2nfS(p A + p B ) = J2 (PA, »"! PB, S\<p) [Xr <8> Xs] , (21) 

and cr^ (erg) is understood to act on % r (% s ) only. 

2.3 Quantum electrodynamics 

As our final example, we turn to quantum electrodynamics (QED). We consider the Coulomb 
gauge for the massless gauge field A which we couple to two different Dirac fields ip A and 
ip B with opposite electric charges ±e. Then the interaction Hamiltonian is 



H 1 = ejd 3 x: [^(x) 7 ^a(x) - $ B (x)-fif> B {x)] Aj(x) : 



e 2 



_ , d3xd3y . [^(x)7Vi(x)-fe(x)7°fe(x)] [Myh°My)-Myh°My)] . 
87T,/ l x ~y| 



(22) 

where A tr denotes the spatially transverse part of the gauge field, the dyamical degrees of 
freedom remaining from A after the Coulomb gauge fixing. They are most simply charac- 
terized by 

A* r (k) = , k ■ A tr (k) = (23) 

for the Fourier coefficients A^ r (k) of A^(x), and can also be projected out of A by applying 
(the negative of) the transverse Kronecker delta <5 tr (k), 

C(k) = ag(k) = ^(k) = o , ag(k) = % - ^ . (24) 

The Bloch- Wilson Hamiltonian for states with one A and one B fermion can be calculated 
in strict analogy to the calculations in Refs. El for the Wick-Cutkosky model and Yukawa 
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theory. Details will be given in a future publication. The result for the matrix elements of 
H B is 

(Pa, r; p B , s\H B \p' A , r'; p' B , s') 

= (V™a + Pa + ym| + p|J (2tt) 3 5( Pa - p^)5 rr .(27r) 3 5(p B - p B )<W 



2E A 2E B 2E A , 2E B , 

PA PB P' A P'g 



X 



1 



(Pa - Pa) 
1 



^2 [uA{pA,r)^°u A (p'A,r')] [mb(pb,s)7°m b (p / b ,s / )] 
1 1 



+ 



2\p A - p' A \ \Eg A + \p A - p' A \ - E$ E$ B + |p fl - p' B \ - El 



x 



[uA(PA,r)ru A (p' A ,r')] I]4 A) (PA - P' A ) 4 A) *(PA - Pa) 



,A=1 



x [u B (p B ,s)Yu B (p' B ,s') 



(2ir) 3 6(p A + p B -p' A - p' B ) 



(25) 



where eW(k) for A = 1,2 are the spatially transverse polarization vectors, so that 

E4 A) (k)4 A) *(k)=^(k). (26) 

A=l 

The Schrodinger equation corresponding to the Bloch- Wilson Hamiltonian (J25)) in the 
center-of-mass system is written in Pauli spinor form as 



m A + p 2 + J m B + p 2 I 0(p 



„2 / P 



E A + m A El 3 + m B fd 3 p> E A +m A E B ,+m B 



2E A 2E B 
p p 



2E$ 2E B , 



1 



+ 



(p - pO 
i 



i + 



p • er A p cr A 



E A + m A E A , + m A 



+ 



1 + 



p ■ a B p ■ er B 



E B + m B Ep, + m B 



2\p-p'\ yE£ + \p-p'\-E$ E B + \p-p'\-E B 

p-a A )a A a\{p' ■ a A )\ t , / (p • (T B )a B a J B (p' ■ tr B ) 



E A + m A E A , + m A 



Eg + m B Ep 3 , + m B 



(p) = E<p(p) , 
(27) 
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with the spinorial wave function 0(p) defined as in Eq. (|2ip. 

3 The Okubo Hamiltonian 

After these concrete examples, let us come back to the general formulae for a moment. 
Our choice of Ub and Hb has the advantages of relative simplicity and the wave function 
interpretation of the ife-eigenstates. However, since Ub is not unitary (this is maybe clearest 
for its inverse, Po\ n ), the Bloch- Wilson Hamiltonian Hb is in general not hermitian, as we 
will see explicitly in the next section for our specific examples. The nonhermiticity of the 
effective Hamiltonian might be a serious drawback in practical applications, although it must 
be mentioned that it has not led to any problems in the numerical calculations of Refs. [21 IS]- 
In any case, a simple possibility is to replace Ub by its unitary part 

U w = U b (U b UbT 1/2 : O ^ Q (28) 

which maps Qo to the same subspace Q as does Ub- Indeed, 

U w = (U B U B )- 1/2 U b = (U B U B ) 1/2 U b x = U w ' . (29) 

The corresponding effective Hamiltonian 

H w = U w l HU w = {U b U b ) 1/2 Hb{U b UbY 1/2 (30) 

is then hermitian, as wished. We will refer to Uw and Hw in the following as the Okubo 
map and the Okubo Hamiltonian to distinguish them from Ub and Hb- They were first 
introduced by Okubo JH] i n a way that does not refer to the generalized Gell-Mann-Low 
theorem or the adiabatic evolution operator. As a consequence, there is no ie-prescription 
for the energy denominators in Okubo's original formulation, and the relation to Feynman 
diagrams is not obvious. 

The Okubo map was first applied to field theory by Wilson in Ref . ^3] > who however did 
not build upon Okubo's work. Later, Gari et al. and Kriiger and Glockle used the Okubo map 
introduced in Ref. [T^j, however without any reference to Wilson's earlier work. To complete 
this historical interlude, it was apparently realized by Wilson's collaborators that the Okubo 
map generalizes the formulae developed by Bloch for degenerate perturbation theory [12j. 
Independently, the present author generalized Bloch's formulation to the Bloch- Wilson map 
and subsequently realized its connection with a generalization of the Gell-Mann-Low theorem 
PP, before becoming aware of the older literature. 

We will now derive the explicit second-order expression for H w analogous to Eq. for 
H B , for the case of particle number changing interactions so that Eq. (JSJ) holds. Equation 
(HJ) implies, due to the unitarity of U e , that 

u B u B = (Pot^r 1 (iw.nr 1 • (3i) 
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(32) 



The expansion (J2J) then leads to 

UhU B = P Q + f dt x I 1 dt 2 e-<^+^P H x (t x )H x (t 2 )P 

J — oo J —oo 

+ f dt x r dt 2 e-<^ + ^P H x (t 2 )H x (t x )P 

J —oo J —oo 

= P + f dt x f dt 2 e-<^ + ^P H x (t x )H x {t 2 )P , 

J —oo J —oo 

to second order in H%. Now, from Eq. (|3L)p we have for H w , again to second order, 
H w = H P - i j dte-^P H x (0)H x (t)P 

J —oo 

+ \ f dh f dt 2 e-<^ + ^Po[H x (t x )H x {t 2 ),H }Po, (33) 

^ J —oo J —oo 

where Eqs. (JHJ) and (|52*j) have been used. The commutator in this expression can be rewritten 



as 



[H x (t x )H x {t 2 ),H ] = [H x {t x ),H ]H x {t 2 )+H x {t x ) [H x {t 2 ),H ] 

= i^-H x (t x )H x (t 2 ) + iH x (t x )^-H x (t 2 ) , (34) 
and integration by parts and the limit e — > finally lead to 

H W = H P - 1 - f dte-^P H l (0)H 1 (t)P + ^- f dt e'^ P H x {t) H x (0)P 
* J —00 J —00 

= ±(H B + H B ), (35) 

i.e., Hw is precisely the hermitian part of Hb (to this order). 

Concerning these calculations, let us remark that there are two formally different repre- 
sentations for the inverse of Uw From Eqs. (J2J), (|2*8)l and ([3*2^1. we have 

fO fO ft! 

U w = P -t / dte-^H x (t)P - dt x dt 2 e- e ^ + ^H x (t x )H x (t 2 )P 
J —00 J —00 J —00 

+ \f dt x T dt 2 e -(N+M)p [H x (t x ), H x (t 2 )} P , (36) 
^ J —00 J —00 

to second order. Then the second-order expression for is different from the formal 
second-order expression for the inverse, the latter turning out to be 

1 r° r° 

U w l = p o+~ dt x / dt 2 e- e(|tl|+|t2|) Po#i(ti)#i(t 2 )Po • (37) 
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However, it is straightforward to verify, by working consequently to second order, that both 
expressions lead to the same result when applied to an element of Q = Ub^o , as they 
must according to the unitarity of Uw, Eq. fP§|) . As a result of the i7-invariance of Q, the 
corresponding expressions for Hw are also formally identical to second order, 

H w = U W HU W = U^HU W , (38) 

hence the present formulation is consistent. The calculation of Hw above implicitly uses the 
expression (pT7j). Note, however, that identities like Eq. (J54*j) are required to establish Eq. 
(they are also required to establish the if-invariance of Q, i.e., the generalized Gell-Mann- 
Low theorem), and that an equality like Eq. (|38|) does not hold for the H - or ii^-part of 
H w individually. 

Equation (|35|) allows to calculate the matrix elements of Hw and the corresponding effec- 
tive Schrodinger equation very easily given the corresponding results for Hb- In particular, 
for all theories considered here, the matrix elements of Hb for the zero- and one-particle 
states are diagonal and real, hence the results for the vacuum energy and the renormalized 
masses including the lowest-order radiative corrections, are identical for Hb and Hw- For 
the Wick-Cutkosky model, the only change in going from Ub to Uw hi Eqs. (jllj) and (jl4j) . 
is the replacement 



2\p A - p' A \ \E£ A + \p A - p' A \ - M Eg B + |p fl - p' B \ - E* 



Pa P s ^ Ps 
1 



+ ,, 2 -7T2 (39) 



2 \(PA~ V' A ) 2 ~ {Ei A - E- a Y (p B - p'bY - {E* B - E* b ) 2 

of the effective potential with its "symmetric part" . In the case of Yukawa theory, one has 
to take into account, in addition, that for the spinor products 

[u A (p' A , r') u A (p A , r)]* = u A (p A , r) u A (p' A , r') (40) 

holds. As a consequence, Eqs. ffTTj) and (|2U|) are unchanged for Hw, except for the replace- 
ment (}3T?j) . The same is true for Eqs. and in Coulomb gauge QED, due to an 
analogous property of the spinor products and the corresponding property of 5 tT [see Eq. 

mi 



4 Expansion around the nonrelativistic limit 

It was emphasized in Refs. j2 El that it is very simple to obtain the nonrelativistic limit 
from the Bloch- Wilson Hamiltonian, by just considering the leading terms for small rel- 
ative momentum. In all the theories considered in the previous section, one obtains in 
this way a nonrelativistic Schrodinger equation with a Coulomb potential, as a consequence 
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of the massless exchanged boson. In the present section, we will show that the next-to- 
leading terms in a systematic expansion in powers of the relative momentum can be used 
to obtain the (lowest-order) fine and hyperfine structure of bound states by application of 
time-independent Rayleigh-Schrodinger perturbation theory. 

The lowest-order fine and hyperfine structure is of course well-known for the case of QED. 
It was first obtained for hydrogenic systems [5] from the Breit equation fBJ. Note that the 
effective Schrodinger equation (|27j) presented in the previous section does not properly apply 
to positronium because there is an additional contribution (the virtual annihilation graph) 
to the Bloch- Wilson Hamiltonian in the case of a fermion-antifermion bound state. Nor does 
it apply to hydrogen because of the anomalous ^-factor of the proton, while g = 2 is implicit 
in Eq. (|27|). It is, however, appropriate for the description of muonium, an antimuon-electron 
bound state. For the other two cases, the Wick-Cutkosky model and Yukawa theory, the 
fine and hyperfine structure have, to our knowledge, not been determined before. 



4.1 Momentum expansion 



The crucial new feature arising in the fine structure for the Wick-Cutkosky model and 
Yukawa theory is the contribution from the retardation in the effective potential. Explicitly, 
the lowest orders of a systematic expansion in powers of the relative momentum are, for Hb, 

p + p' p — p' 



|p-p'| E£+\p-p>\-E$ 



IP-P' 



2m 



A |P-P' 



p + p' p — p' 

2m A |p - p'| 



+ O [{p/m A f] 



(41) 



where p = |p|. The leading term corresponds to the Coulomb potential, while the first- 
order, antihermitian correction term and the second-order hermitian term will turn out to 
contribute to the same order of the fine structure. On the other hand, for H\y only the 
hermitian part contributes, explicitly 



1 



1 



(P-Pf -(^"^) (P-P' 



p + p' p — p' 
2m A |p-p'| 



+ O [(p/m A ) 4 ] 



(42) 



Since the antihermitian term in Eq. (j41|) does contribute to the lowest-order fine structure, 
the effective Hamiltonians Hb and Hy/ lead to different results to this order, for the Wick- 
Cutkosky model and Yukawa theory. This fact and its possible implications will be discussed 
further towards the end of this section and in the next section. 

We will now present the complete results for the expansion in powers of relative mo- 
mentum to next-to-leading orders, for the different effective Hamiltonians discussed in the 
previous section. For the Wick-Cutkosky model and the Hamiltonian Hb, the expansion of 
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the effective Schrodinger equation (fTlj) to next-to-leading order gives 



2m r 



im. 
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b / J 
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d 3 p' 
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p'Y 
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■IP - P'l 
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p - p') 8m| (p - p')' 



P 



(43a) 
(436) 
(43c) 

m A - m B )4>{p) , (43d) 



where we have introduced the reduced mass 



m A + m B 

Analogously, for the Yukawa theory and the Hamiltonian H B , Eq. (f2T)|) leads to 



2m T . 
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B/ J 



0(P) 
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d 3 p' 

4m r |p — p 
(p 2 — p' 
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+ 



^2 _|_ p/2 



+ 



(p2 _ p /2)2 



8ra^ (p — p') 8m| (p — p' 
(P ■ 0"a)(p' • <Ta) (p • cr B )(p' ■ <x B ) 



4m^ 



4m B 



(p') = (E -m A - m B )(j)(p) 



(44) 

(45a) 
(456) 
(45c) 
(45d) 
(45e) 



to next-to-leading order. Note the differences between the two effective Schrodinger equa- 
tions, a relative factor of 1/2 between Eqs. (|436J) and (J456)) . and the appearance in Eq. (|45ejl 
of the Pauli matrices acting on the Pauli spinors of 0(p') [see Eq. (|21j) ]. both originating 
from the product of Dirac spinors in Eq. (fTTj) . 

Turning now to QED in the Coulomb gauge, we can see directly from the effective 
Schrodinger equation (|27|) why there is no contribution from the retardation to the lowest- 
order fine structure: the instantaneous Coulomb interaction of the charge densities has no 
retardation, and the interaction with retardation transmitted by the spatially transverse pho- 
tons carries additional powers of momentum, so that retardation only contributes to higher 
orders. As a consequence, there is no antihermitian term in the expansion to next-to-leading 
order, and H B and Hyy lead to the same fine and hyperfine structure. 

In order to obtain a compact explicit expression for the next-to-leading terms in the 
expansion, some algebra involving the Pauli matrices is required. In particular, the identity 



(p • cr A )a A + (T A (p ■ (T A ) = p + p + i er A x (p - p 



(46) 
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is helpful in intermediate steps. The final result can be written as 
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(p-o-a)(p / -o-a) (p • <Tb)(p' ' o-By 



IrriAvriB 
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AmATTlB 



(p - p 



/\2 



Ivhaviib 
(p 2 -p' 2 ) 2 

(p - pO 2 



(p') = (E-m A - m B )0(p) . 



(47a) 
(476) 
(47c) 
(47d) 
(47e) 
(47/) 



Expressions ()47a|) , (j476|) , and ()47c|) are equal to ()45a|) , ()45b|) , and ()45e|) , except for the oppo- 
site sign of ()47cj) compared to ()45ej) . However, as discussed above, there are no contributions 
in the QED case analogous to (|45cj) and (|45<i|) from retardation, while in Yukawa theory there 
are no terms analogous to (|47rfjl . f)47ej) . and J47/D - It is interesting to remark that the latter 
terms tend to zero in the one-body limit m B — > oo, so that the lowest-order fine structure 
calculation in the one-body limit is much simpler for the (well-known) QED case than for 
Yukawa theory. 

Now, if we consider the effective Hamiltonian Hy/ instead of Hb, according to the results 
of the previous section and Eq. (|42j) . the only change in the effective Schrodinger equations 
is the absence of the antihermitian terms (|43cj) for the Wick-Cutkosky model and (|45cj) 
for Yukawa theory, while the effective Schrodinger equation for QED is left unchanged as 
mentioned before. 



4.2 Perturbation theory for the hermitian terms 

As far as the absolute order of the correction terms is concerned, we note the well-known 
fact that p/m r = 0{a) in the nonrelativistic limit (from the expectation value of p 2 ), where 
a is the fine structure constant. Although we will use the symbol a indistinctly, it is defined 
differently in the three theories we consider: we take 

a = 9 - (48) 

IQ-KniArriB 



in the Wick-Cutkosky model [compare with Eq. (J43 



a = £ (49) 
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in Yukawa theory, and, as usual, 



e 2 



a = — (50) 

in QED. Note that p/m r = 0(a) is consistent with the nonrelativistic result for the binding 
energy 

E — rriA — rriB = 0(m r a 2 ) (51) 

by naive (small-p) power counting when we count powers of p and p' indistinctly and take 
the integration measure d 3 p' into consideration. This kind of momentum power counting 
which we will refer to as "IR power counting" , will be used extensively in Appendix [Bl 
By the same IR power counting, the perturbative corrections contribute to the order m r a 4 
[the contributions from Eqs. (|47(ij) . ()47e|) . and Q47/D are suppressed by a factor m^/m^ 
for vtla <C tub], except for the antihermitian terms ()43cj) and ()45cj) which are formally of 
0(m r a 3 ). However, as mentioned before, the contributions of the antihermitian terms to this 
order turn out to vanish and they rather contribute to 0(m r a A ) like the hermitian terms. 

For the rest of this section, we will discuss the application of Rayleigh-Schrodinger per- 
turbation theory in order to obtain analytical results for the fine and hyperfine structure in 
the different theories. In view of future applications with higher-order perturbative contribu- 
tions, we will perform all calculations directly in momentum space, and not in position space 
like in the traditional QED fine structure calculations. For all hermitian contributions, we 
apply Rayleigh-Schrodinger perturbation theory to first order which is quite straightforward. 
The nonrelativistic wave functions in momentum space needed for the calculations are listed 
in Appendix [X] for completeness. The matrix elements in these angular momentum eigen- 
states are determined by decomposing the perturbations H pert in partial waves and applying 
the well-known spherical harmonics addition theorem, 

00 4tt - 
H peit (p,p', cos 9) = J2^ 1 -a l (p,p') Y lm (p)YL(p) , (52) 

1=0 m=-l 

where p = p/p, 9 is the angle between p and p', and 

, 21 + 1 f 1 

a>i(p,p') = — ^ — J dcos9Pi(cos9)H pcrt (p,p',cos9). (53) 

The explicit expressions for the relevant coefficient functions ai(p,p') are easily determined. 
For the perturbations that involve the Pauli matrices, the formulae developed in Ref. jS] 
for the application of the helicity operators (p • er) to total angular momentum eigenstates 
(including the spins of the two fermionic constituents) are used. For convenience, they are 
also reproduced in Appendix El The remaining integrals over p and p' are elementary, with 
the exception of the integrals 

dp 7 — 5—5^ — 9N „ t i In ( , P + ) , (54) 
(rrrcr + n z p z ) n+L \\P — P\J 
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where n — 1, 2, 3, . . . and < k < n. All these integrals can be obtained by differentiation 
with respect to A 2 and trivial algebraic manipulations of the fraction from 

dp tj-^ — o In ( r~~~^7i ) = 71 arctan(j// A) . (55) 



o \ 2 +p 2 \\p-p'\; 
Equation (|B3j) is derived in Appendix |X] with the help of complex contour integration. 

4.3 Perturbation theory for the antihermitian term 

Let us now develop the analogue of Rayleigh-Schrodinger perturbation theory for antihermi- 
tian perturbations as they appear in the Wick-Cutkosky model and Yukawa theory. Consider 
the general case of a Hamiltonian of the form H = Hq + H a with Ho hermitian and H a an- 
tihermitian. The eigenstates and eigenvalues of H are supposed to be known, 

Ho\€ ) ) = E^). (56) 

Let us assume for simplicity that the eigenvalues (or at least the one considered) are not 
degenerate. This is not true in our case, however, the submatrices of H a in the degener- 
ate subspaces are diagonal and the formulae developed in the following apply despite the 
degeneracy. 

We expand the eigenstates and eigenvalues, as usual, around the ones of H Q , 

I0n) = \€ ] ) + + 

En = E%> + Ej? + E^ + .... (57) 

Since the full Hamiltonian is not hermitian, we cannot a priori assume that the eigenvalues 
are real nor that the eigenstates are orthogonal, and the same is hence true for the correction 
terms in the expansions above. However, we will insist on the normalization of the eigenstates 
and also adopt the usual phase convention 

(*«|*«>) real. (58) 

For i = 1, the above implies that 



Uo)U(i) 



) = . (59) 



Inserting the expansions into the Schrodinger equation for H, one has the infinite tower 
of equations (beginning with (JHHJ) ) 

Ho\€ } ) + H a \€ ] ) = E n °W } ) + E?m , (60) 

H \<f> { : ] ) + H tt \W) = E$>\4>M) + EjpW) + E*>\<I>M) , (61) 

etc. Projecting with (0^ o> | on the first of these equations leads to the analogue of the well- 
known result in hermitian perturbation theory, 

= (W\Ba\W) = -mH a \W , (62) 
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where we have explicitly used the antihermiticity of H a in the second equality. 

In our case, H a in Eqs. (|43c|) and (J45c|) conserves spin and orbital angular momentum sep- 
arately, hence the angular and spin dependence of the wave functions 4> l ° } is unchanged under 
H a , and the diagonal matrix elements reduce to integrals over the moduli of momenta. Since 
the operator H a and the radial (zero-order) wave functions are real, Eq. (}6*2j) implies that 
E^ = 0, i.e., the antihermitian perturbations give no contribution in first-order perturbation 
theory. However, H a is of order p/m relative to the leading term in the nonrelativistic limit, 
while all the hermitian perturbations are of order (p/m) 2 , hence the contributions of H a in 
second-order perturbation theory are potentially of the same order in a as the contributions 
of the hermitian perturbations in first-order perturbation theory. 

To determine the second-order contributions of H a , we use Eq. (JoTij) again, but now project 
with (0^1, m ^ n, to find an expression for ((fi^l^) , and using the completeness relation 
and the phase convention (J59|) one finds again for the first corrections to the eigenstates the 
analogue of the result in hermitian perturbation theory, 



W) = E IC VS^'ff > (63) 

mpn 

where the sum runs over all zero-order eigenstates \4>^) with m ^ n. 
Finally, project with (0^ o) | on Eq. (JBTj) to find 

E*> = (W\H a \4>*>) , (64) 

from where it follows by use of Eq. that 

mpn mpn 

The sign appearing in the last step stems, of course, from the antihermiticity of H a . It 
implies, in particular, for the correction to the ground state energy, £q 0) < E^ for all 
m 7^ 0, that Eq > 0, contrary to the case of a hermitian perturbation. Note that E^ } 
turns out to be formally of 0(m r a 4 ), as anticipated, if we take (^\H a \^) = 0(m r a 3 ) and 
= 0(m r a 2 ). 

Now the sum in Eq. (|65|) is not always simple to evaluate analytically. We will procede 
here in analogy with the method of Ref. |T7j. The idea is to use Eq. (J6U)) directly to determine 
\<f>$) rather than employing the expansion (|63p. The correction E^ is then found from Eq. 
(jMJ). In our case, Eq. (j6TJj) for \<p^) takes a particularly simple form due to the fact that 
E^ = 0, namely, 

H W) = E^) - H a \<f>W} , ( 66 ) 

which is the Schrodinger equation for the eigenvalue E^ (with solution |0,^ o) )) with an ad- 
ditional inhomogeneous term. What follows is the only part of the calculation that we have 
performed in position space, for reasons to become clear shortly. It is not difficult, of course, 
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to transform every step of the argument to momentum space. We will also present a differ- 
ent approach to the calculation of in Appendix [B] which proceeds entirely in momentum 
space. 

To proceed with the calculation in position space, we then need the explicit expression 
for H a in position space. To this end, start with 

' 1 r d 3 re- lip - p ' ) r \ . (67) 



|p — p'| 2tt 2 J r 2 
Taking the derivate of Eq. (|67|) with respect to p and multiplying with (p + p') gives 
p 2 — p' 2 1 f , r / d A r j , d 



p — p'| 3 2tt 2 J r 2 \ dr dr 



d 3 r--( e ~ ipr ^e tp r - e ip r ^ e - ipr , (68 



so that 



(69) 



Integrating over p, p', and r' finally leads to 



r> /2 8 1 

Ha^(r)-^-r + (70) 



2nm r \r dr r 2 / 
For the solution of Eq. (jHEJ), we now make the ansatz 

€L(r) = Xmir'Wjr) , (71) 

where 

CUM = ^Y lm (i) (72) 



r 

is a solution of Eq. (jSEJ), and still has to be to be multiplied with the appropriate Pauli 
spinors to describe the spin orientations of fermions A and B in the case of Yukawa theory. 
When we insert this ansatz into Eq. (JBljj) and use the Schrodinger equation (Jo7)j) for <^°/m( r )' 
the equation 



for Xni{ r ) results. The terms on the left-hand-side of Eq. (|7Sj) stem from the second r- 
derivative in the kinetic term of Hq. 

Equation (fTHJ) is identically fulfilled, i.e., independently of the function u °i(r), for 

d-Xni & (74) 

dr irr 
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Being a linear differential equation for dx/dr, Eq. (|73|) has a one- dimensional (affine) vec- 
tor space of solutions. We will hence have a look at the solutions of the corresponding 
homogeneous equation. Neglecting powers of r, u { °l(r) falls to zero for large r like 

uS'CO-e-*/™*, (75) 

do being the Bohr radius. It follows that the nontrivial solutions of the homogeneous equation 
behave like 

<bU _ e 2r/na ^ 

dr 

for large r, and the corresponding functions 0„/ m (r) in Eq. (fTTj) are not acceptable as bound 
state solutions. As a result, we have to take the trivial solution of the homogeneous equation, 
and Eq. (jTlj) is in fact the physical solution of Eq. (fT3j) . 

The simplicity of the ansatz (j71)) which leads to the solution (|74jl . is the reason for our 
present use of position space. In momentum space, Eq. (fTTj) corresponds to a convolution 
which makes the latter formulation a little less transparent in the solution of Eq. ()66j). 

From the integration of Eq. (fTlj) we obtain 

ot 2r 

Xni{r) = -\n + C n i, (77) 

7r na 

with a convenient normalization of the argument of the logarithm. The integration constant 
C n i is fixed by the phase convention (foTIj) which is where the (n, /)-dependence of Xni( r ) enters 
[compare with Eq. (|74~j) ]. The constant actually plays no role at all in the calculation of the 
perturbed energy with the help of Eq. (jMjl . given that ((p^lHal^) — as discussed after 
Eq. ((021)- 

The constant C n \ is important, however, when one is interested in the wave functions 
themselves. The integrals corresponding to Eq. (foTJj) can be evaluated with the help of 
integral tables, for example, Ref. [TH] - The explicit results for the lowest-lying states are 

a ( 3 



7T V 2 

a / 9 

= — I lE 7 



7T \ 4 

a ( 25\ 

C 2 i = - \1E - , (76 



71 v 12 y 

with the Euler-Mascheroni constant 7^. It is tempting to speculate that the logarithms of r 
in the corrections to the wavefunctions, when summed up to all orders of the perturbative 
expansion, combine into an a-dependent power of r. This gives a hint to a possible improve- 
ment in the choice of the basis functions used for the numerical solution of Eqs. (|14j) and 

(HD. 

From Eqs. (JMj). (fTTj) . and (|77jl we can now calulate the corrections E n 2) to the energy, 
again with the help of integral tables. However, there is a different, interesting, and somewhat 
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simpler way to calculate these corrections. To get there, we manipulate the expression 
for in the following way: 

£W _ « r drr^ {T) ( 2d I ^^WW 



2iim r J r \r dr ' r 2 

" ^ drr 2 ^^ (2,h ^ ""' ir) 



27im r J r \r dr 



, g r dTr 2 X{r)<\r) (2 d | 1^ <»(r) 



2iim r J r \r dr ' r 2 



2Txm T ^ n 1 nr 2 n 2-Km r J r \r dr r 2 / 

where we have used Eq. (|73j) in the last step. It follows that 

which has the form of a first-order correction for the hermitian perturbation 

c? 1 



27i 2 m r r 2 

or, equivalently, 

a 2 1 , . 

— i r 82 

m r |p — p | 

in momentum space. Expression (|50jl can then be evaluated by the standard methods 
desribed in the previous subsection. 

Let us emphasize again that the contribution (JHUj) is due to the retardation of the interac- 
tion in the Wick-Cutkosky model and Yukawa theory, and that it is absent in the hermitian 
Okubo Hamiltonian. This situation is reminiscent of a discussion in the older literature on 
the effective one-boson exchange (OBE) description of the nucleon-nucleon interaction: the 
Gross equation develops a repulsive contribution to the interaction potential in a (not com- 
pletely systematic) expansion around the nonrelativistic limit, while such a contribution is 
absent in the Blankenbecler-Sugar-Logunov-Tavkhelidze (BSLT) equation [19] . Incidentally, 
the Gross equation includes retardation in its kernel while the BSLT equation does not (to 
lowest order). In Ref. jT9], the repulsive contribution to the potential was associated with 
the repulsive core of the nucleon-nucleon interaction. Since the exchange of a scalar plays 
an important role in the effective OBE description giving the dominant attractive contri- 
bution to the intermediate-range potential, it is actually possible that the contribution ([ST]) 
is relevant to physics in the context of the OBE potential. Another interesting remark in 
this respect is that the full Gross equation is plagued by singularities which complicate its 
numerical solution [20], while there is no such problem with Eq. (|2U|) 3j. 
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How is it possible that the two effective Hamiltonians Hb and Hw lead to different 
predictions about the lowest-order fine and hyperfine structure? Both effective Hamiltonians 
have been calculated to order H\ or g 2 , and the next terms in the perturbative expansion 
are of the order g A . It is possible that these latter terms contribute to order m r a 4 , i.e., to 
the lowest-order fine and hyperfine structure, and these contributions need not be the same 
in the two cases. Summing all contributions to order m r a 4 , from the order-g 2 and order-g 4 
terms of the effective Hamiltonian, however, must lead to the same (the complete) result 
for the lowest-order fine and hyperfine structure. It is then also clear that we cannot be 
sure that the terms of order a 4 obtained in the present paper represent the complete fine or 
hyperfine structure in any of the two cases. In fact, the difference between the results for 
Hb and Hy/ shows that at least one of them (or both) must be incomplete. From the results 
of Appendix |B1 it appears plausible that the contribution (|HT|) is obtained from the order-g 4 
term of the effective Hamiltonian Hw, so that at least in the latter case the complete fine and 
hyperfine structure probably cannot be obtained from the g 2 -term alone. On the other hand, 
the g 2 -term is sufficient to generate the complete fine and hyperfine structure in the case of 
QED (in the Coulomb gauge), for both Hb and Hw, as we will show in the next section. In 
particular, then, there really are cases where the g 2 -term is sufficient for our purpose. The 
situation remains unclear in the Wick-Cutkosky model and Yukawa theory until we have 
determined all terms of the order g 4 in the effective Hamiltonian, although it might well be 
that the g 2 -term in Hb (and thus not in Hw) generate the lowest-order fine and hyperfine 
structure completely. In fact, it may be a reasonable criterium for a "good" relativistic 
bound state equation that the lowest-order approximation produce all the lowest-order fine 
and hyperfine structure. 

5 Results and discussion 

We will now give explicit results for the different contributions to the fine and hyperfine 
stucture in the three theories considered, for the states with principal quantum numbers 
n = 1,2. These results will be discussed with respect to consistency, comparison with the 
numerical solutions of the full effective Schrodinger equations and with other relativistic 
bound state equations. In the case of QED, the well-known results will be reproduced 
completely, while the results for the lowest-order fine and hyperfine structure in the Wick- 
Cutkosky model and Yukawa theory have, to the best of our knowledge, not been obtained 
before. 

5.1 Wick-Cutkosky model 

By applying the methods detailed in the previous section (and in the appendices) to the 
lowest-lying states, we arrive at the explicit results for the relativistic energy corrections 
AE(nL) (L = S,P,D,.. .) cited below. The different contributions specified follow the or- 
der in the expanded effective Hamiltonian Hb, i.e., from left to right we have the relativistic 
correction to the kinetic energy in Eq. (|43ajl . the correction term in Eq. (j436|) . the antihermi- 
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tian term ()43c|) . and the term on the left-hand side of Eq. (|43rf|) . The corresponding results 
for Hw are obtained by just leaving out the contribution from the antihermitian term (third 
term in every result). The explicit expressions are 
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(83) 



The most important general features of these results are the following: the sum of all cor- 
rections is positive for any of these levels, and for any mass ratio m^/wi^. After taking into 
account these corrections, the energy of the 2S level lies above the 2P level (these levels 
are, of course, degenerate in the nonrelativistic limit). Both of these features are due to 
the dominant correction term from Eq. (j436|) . This term arises from the expansion of the 
inverse square roots of kinetic energies in Eqs. (jllj) or (J14)) . which are characteristic of the 
nonlocalizability of point particles in relativistic quantum theories. The overall contribution 
of retardation, corrections (|43c|) and (|43dj) . is negative, hence the retardation of the potential 
has an attractive effect (see, however, the discussion of the antihermitian term at the end 
of Subsection 14. 3|) . Incidentally, the relation between the energy corrections from these two 
terms is the same for all energy levels. 

The most important property of all, however, is probably the fact that the corrections 
to the energy are of order m r a 4 , as one might have expected from the well-known fine 
and hyperfine structure of hydrogen. A look at Fig. ^ which we have reproduced here, 
for convenience, from Ref. [2], shows that this feature is nicely reproduced by the results 
from a numerical solution of the full effective Schrodinger equation (fT4"|) [2]. Observe that 
the energy eigenvalues presented in Fig. are normalized to the nonrelativistic inonization 
energy m r a 2 /2, hence the order m r a A of the relativistic corrections corresponds to a quadratic 
behavior of the curve near a = 0. 

It is a remarkable fact that the numerical curves for all other relativistic bound state 
equations presented in Fig. ^ do not show a quadratic behavior for small a (disregarding the 
results from the Haag expansion which does not have a proper nonrelativistic limit), which 
points to the existence of an order- (m r a 3 ) term in their expansion around the nonrelativistic 
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Figure 1: The spectrum of binding energies £?& = E—rriA—iriB for 5— states in the equal-mass 
case, compared to the ground state energies of the Wick-Cutkosky model [TU], the Hamil- 
tonian eigenvalue equation in a Fock space truncation the Regge theory predictions 
[22] j the light-front calculation [23] , and the Haag expansion results (21j in their domain of 
validity. 

limit. In fact, it is known that there appears a "curious" term of order m r a 3 In a in the 
ladder approximation to the Bethe-Salpeter equation in such an expansion [11] . A simple 
power-counting argument cannot exclude the existence of contributions to the order m r a 3 
(or m r a 3 In a) from the terms of order Hf in the expansion of Hb derived from the Dyson 
series (J2J) which are not considered in the present paper. However, from our experience 
with QED (e.g., the hydrogen atom and positronium) including terms of the order Hf, such 
contributions are certainly not expected. Note that a contribution of the next-higher order 
m T a A from the order- Hf terms is necessary for the consistence of the results for Hb and H\y 
presented here which differ by a contribution of the order m r a 4 , as we have seen (compare 
with the discussion at the end of Subsection 14. 3|) . The definite solution of these issues, most 
importantly the question of whether the results presented here for the fine structure (from 
Hb) are complete, has to await the (complicated) evaluation of the O(Hf) terms by the 
present or a different method. 

To give a better idea of how the numerical results relate to the perturbative expansion 
around the nonrelativistic limit, we have plotted both in Fig. Elfor the ground state. The 
approximation by the perturbative expression is good only way below a = 0.1. We have also 
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Figure 2: The numerical results (diamonds) for the ground state binding energy, normalized 
as in Fig. [TJ compared with the perturbative results (JH3|) (solid line). We have also plotted 
the perturbative results without the contribution from the antihermitian term (dashed line). 

plotted the curve corresponding to all hermitian terms, i.e., omitting the contribution from 
Eq. ()43c|) in Eq. (JHHj) . Note that the latter curve lies closer to the numerical data points than 
the curve for the complete expressions, for intermediate values of a. 

5.2 Yukawa theory 

For Yukawa theory, the calculation is technically, but not essentially, more complicated than 
for the Wick-Cutkosky model, due to the spin of the fermions. The necessary techniques 
for the evaluation of the relativistic corrections are presented in Appendix [XJ We will 
label the states by the standard spectroscopic notation n 2S+1 Lj. The explicit results for 
principal qunatum numbers n — 1,2 are, again term by term in the order they appear in the 
momentum expansion of the effective Hamiltonian, Eqs. ()45a|) - (j45e|) . 
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The results for the states 2 3 Pi and 2 x Pi bear primes because these states (which are de- 
generate in the nonrelativistic limit) mix through the equal off-diagonal matrix elements 



x/2 
96 



m, 



"A m B, 

To obtain the energy corrections properly, one hence has to diagonalize the corresponding 
2x2 matrix with diagonal elements AE'(2 3 Pi) and AE'(2 1 Pi). The explicit expressions for 
the eigenvalues are not too illuminating. 

Compared to the results for the Wick-Cutkosky model, in the Yukawa case the contri- 
butions from the "nonlocal" term (|436|) are reduced to half their value by a contribution 
from the normalization of the Dirac spinors in (j45fe|) . and a positive contribution from the 
spin structure enters. The fact that the "nonlocal" contribution is reduced leads to a sign 
change for the relativistic corrections of the IS* states depending on the mass ratio: in the 
case of equal masses, all corrections are positive, while the corrections for l^So and l 3 5*i turn 
negative in the one-body limit mg-* oo. 

As for the "mixing" states 2 3 Pi and 2 1 Pi, in the case of equal masses the off-diagonal 
elements (18 5 j) vanish, and AP"(2 3 Pi) and AE'(2 1 Pi) directly give the energy corrections of 
the states. Observe the equal spacing between the energy corrections of the four 2P states 
in this case which is nicely reproduced by the numerical results up to intermediate coupling 
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constants. In the one-body limit, on the other hand, the states mix strongly, and one of the 
eigenstates becomes degenerate with the state 2 3 P , the other with 2 3 P 2 . These degeneracies 
and the ones of l 1 ^ with l 3 ^ and of 2 1 S , with 2 3 S' 1 , are exact in the one-body limit, for 
any value of the coupling constant. The reason for these twofold degeneracies is that the 
spin of fermion B decouples from the dynamics in the limit -* oo. In fact, the states 
2 3 P and the mixture of 2 3 Pi and 2 1 P i degenerate with the former, are eigenstates of the 
total angular momentum (relative orbital angular momentum and spin) of fermion A with 
eigenvalue ]a = 1/2, while state 2 3 P 2 and the other linear combination of 2 3 P 1 and 2 1 P 1 are 
2a eigenstates with eigenvalue 3/2. Note that the ja = 1/2 state lies higher in energy than 
the ]a = 3/2 state. 

Close to the one-body limit, we can check for hyperfine structure by expanding the 
expressions for the energy corrections of the 2P states in Eq. and the off-diagonal 
elements (J85)) in powers of m^/m^. Obviously, there is no hyperfme splitting for the IS" 
and the 25' states to any order in m^/^B (and to order m r a 4 ). After diagonalizing the 
submatrix for the states 2 3 Pi and 2 l P\, one finds no hyperfine splitting to the order m^/f^B, 
while there does appear such a splitting to the next order, (mi/mg) 2 . These results on the 
hyperfine structure have been anticipated in Ref. [3]. 

The perturbative results (|84j) show all the qualitative features that we have found in the 
numerical solutions of the effective Schrodinger equation (I2U|) in Ref. j^J. In quantitative 
terms, the agreement between numerical and perturbative analytical results is similar to the 
case of the Wick-Cutkosky model (see Fig. EJ) for the S states, and somewhat better for 
the P states. Just like in Fig. the perturbative curve without the contribution of the 
antihermitian term lies closer to the numerical results for intermediate coupling constants 
than the full perturbative curve. It appears that the antihermitian term is only important 
for very small a and does practically not contribute for larger values of a. The most clear- 
cut case is the one of AE(2 1 Sq) = AP(2 3 Si) in the one-body limit where the complete 
perturbative correction is positive, while its hermitian part is negative. In this case, the 
numerical values for very small a increase with the coupling constant as expected from the 
perturbative results, but then start to decrease and become negative roughly at a = 0.35, 
loosely following the "hermitian" curve for larger a. 

The peculiar role that is played by the contributions from the antihermitian term might 
lead one to the believe that these contributions are spurious and may cancel with other 
contributions from terms of the order Hf in the expansion of the effective Hamiltonian. In 
this sense, the Okubo Hamiltonian Hy/ would be a better choice for an effective Hamiltonian, 
because it simply does not contain the antihermitian term. However, the results of Appendix 
iBl suggest that, quite to the contrary, the terms of the order Hf in the expansion of Hw 
generate the hermitian equivalent of the term in question. 
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5.3 Quantum electrodynamics 



Finally, we will present our results for QED in the Coulomb gauge, i.e., the expanded effective 
Hamiltonian of Eqs. (|47aJl-(4TJ). The explicit expressions are [22], term by term, 
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As in the case of Yukawa theory, the states 2 3 Pi and 2 x Pi mix through the equal off-diagonal 
matrix elements 



\/2 / J 1_ 

96 \ m /i m ? 
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M " b Bj 

When we compare these results with Eq. ()84j) for Yukawa theory, the differences are the 
opposite sign of the "spin-orbit" term [Eq. (|45e|) vs. Eq. (|47c|) ]. the absence of retardation, 
and, of course, the terms ()47rf[)— ( 47j I which introduce the hyperfine structure. The difference 



in signs of the spin-orbit terms changes the qualitative features of the spectrum drastically: 
nearly all energy corrections are negative, only AP(l 3 Si) becomes positive for mass ratios 
itla/^b close to one. The level ordering of the 2P states is just opposite to the Yukawa case. 
In the one-body limit, the characteristic degeneracies appear, but in addition the 2S states 
become degenerate with the 2P states with ja = 1/2 (compare with the discussion of the 
one-body limit in the previous subsection). This is, of course, just the famous degeneracy in 
orbital angular momentum characteristic of the Coulomb potential (note that this degeneracy 
is broken in Yukawa theory by the retardation terms). Also, the 2P states with ja = 1/2 
have a lower energy in QED than the ]a = 3/2 states. 

The hyperfine structure terms split the energies of the IS* and the 2S states. Close to the 
one-body limit, we can again expand the expressions (J86|) in powers of m^/mg to obtain the 
hyperfine splittings (after diagonalizing the 2x2 matrix for the states 2 3 Pi and 2^1), which 
now appear to the order mi/rag (and m r a 4 ). Our results (|56^l coincide completely with 
the ones in Ref. jS] close to the one-body limit and for the case of equal masses [23] (in the 
latter case, they coincide with the positronium results omitting the contributions from virtual 
annihilation there). In fact, one can show that the expanded effective Schrodinger equation in 
Eqs. (I47a|) -( |477| ) is identical to the Pauli approximation of the Breit equation (and in the one- 
body limit to the Pauli approximation of the Dirac equation) [5 . Depending on the form in 
which the latter equations are written, the identification of terms may require some algebraic 
labor. In particular, the Pauli approximation of the Breit and Dirac equations are sometimes 
given in a not manifestly hermitian form. In any case, the identity of the equations implies 
the identity of the lowest-order relativistic corrections, for any state and any mass ratio. 
We conclude that our approach reproduces the lowest-order fine and hyperfine structure 
completely (at least) in the case of QED, from the effective Hamiltonian (Hb or H\y) to 
order Hf, and it does so in a very transparent and economic way. 

To sum up, we have argued in this paper that we can obtain the lowest-order fine and 
hyperfine structure in nearly nonrelativistic bound systems in a very straightforward way 
from the application of the generalized Gell-Mann-Low theorem. We have verified this 
claim for bound states in QED, and we have presented, for the first time to our knowledge, 
results for the fine and hyperfine structure of bound states in the Wick-Cutkosky model 
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and Yukawa theory. Although the results are physically appealing, we have as yet no proof 
that our calculation of the fine and hyperfine structure in the latter theories is complete. 
We have emphasized the contributions from the retardation in the effective potential for the 
Wick-Cutkosky model and Yukawa theory, and discussed the role of a peculiar antihermitian 
term that appears in this context, in particular with respect to two possible choices Hb and 
H w for the effective Hamiltonian. 
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A A collection of useful formulae 

We will present in this appendix the formulae that are used in the actual analytic calculations 
of the perturbative corrections to the Coulomb energy eigenvalues. Since we have done all of 
the calculations directly in momentum space (as a warm-up for more complex calculations in 
the future), we need, first of all, the Coulomb wave functions <fi n i m (p) in momentum space, 
here for the principal quantum numbers n = 1,2: 

167r(m r a) 5 / 2 

0ioo (p) = ~ — — 2 ^oo(p) 

[{m r Q>y + p z \ 



'200 (P) = ^3 loo(p) 

(m r a) 9 
V 4 +P 

47m/- (m r a) 7/2 p 

WP) = „ V 6 r 2 73 *Wp) • ( 88 ) 

,m r a) 9 

- — -—- h p 



The wave functions are normalized to 

d 3 p 



(2nf 



|0nta(p)| 2 = l. (89) 



In the cases of Yukawa theory and QED, the expressions above still have to be multiplied 
with the appropriate Pauli spinors for the spin orientation of the fermions. It is convenient in 
the latter cases to use total angular momentum eigenstates. In order to make contact to the 
usual spectroscopy, we first couple the two spins 1/2 to a total spin S — 0, 1, and then couple 
the total spin with the relative orbital angular momentum I to the total angular momentum 
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Xi-i 



J, M. The explicit results for the total angular momentum eigenstates 2S+1 3^i/(p) are [Sj, 
in terms of the well-known eigenstates Xs,m s °f total spin: 

= Y JM (p) xoo , 
1 

" y/2J(2J-l) 

+ y/2(J-M)(J + M)Yj_ 1<M {p) xio 

+ y/(j + M-i){j + M) y>_i^_i(p) XI 1 
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+ v /(j _ M + 1)( J - M + 2) r J+liM -i(p) Xn 
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For the calculation of the perturbative corrections in these theories, one needs to apply 
the operators (p • <r a,b) to these total angular momentum eigenstates. The corresponding 
formulae were also worked out in Ref. They read 

(p ■ °a) yUi>) = v ^2JTi 3 ^ j -m/(p) - {yTV\ ^."(p) 



(p-o-a) ^-i,m(p) = P 

(p ■ <X A ) = - 

(p ■ cta) 3 ^i + i iM (p) = - 
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and 
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(p ■ *b) 3 yj +1>M (P) = p [ l y J JM {i>) - V^pi %,(p) ] ■ (92) 

In the special case J = 0, the states 3 yj_i m(p) an d 3 3^jm(p) do not exist, and on the right- 
hand sides for the application of one of the helicity operators to ^jm(p) an d 3 3^j+ia/(p)? 
only one term remains. For the application of (a a • &b), we use the well-known identity 

(<TA ■ <tb) 2S+l yL(p) = [2S(S + 1) - 3] 25+1 y/ M (p) . (93) 

After the application of these formulae, the angular integrations can be performed with 
the help of the partial wave decomposition and the spherical harmonics addition theorem, 
Eq. (|52|) . To this end, the coefficient functions have to be calculated from Eq. (p)3*|) which is 
elementary in all relevant cases. 

Finally, the integrations over p and p' have to be performed. All appearing integrals are 
again elementary, the only exception being integrals of the type 



(X 2 +p 2 ) n+1 \\p-p 1 

with n — 1,2,3,... and < k < n. All these integrals can be obtained by differentiation 
with respect to A 2 and trivial algebraic manipulations of the fraction from 

ip^n(l±^), (95) 



A 2 + p 2 \ \p — p' 
which we will now evaluate. To this end, consider the analytic function 

F{p) = ^L_ t jE±£) (96) 

X 2 + p l \P — PJ 

which is even under p — > —p and its real part coincides with the integrand of /, so that 

1 r°° 

I = - Re / dpF(p) . (97) 

2 J -oo 

As far as its analytic structure is concerned, F(p) has simple poles at p — ±zA and a cut on 
the real axis from —p' to p'. We can hence evaluate I by integrating along the real axis and 
(e.g.) the upper rim of the cut and closing the contour through the upper infinite semicircle. 
Then only the pole at iX is picked up (suppose A > 0), and the result is 



7r arctan(p'/A) . (98) 
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B A transformation of Foldy-Wouthuysen type 



In this appendix we will devise an alternative method for the calculation of the contributions 
to the energy corrections that arise from the antihermitian term ()43c|) or ()45c|) in the Wick- 
Cutkosky model or Yukawa theory, respectively. Unlike the "direct" calculation in Section 
14 .3| we will construct a similarity transformation that converts Hb to a hermitian operator, 
to the perturbative order required. Since the eigenvalues are unchanged under the similarity 
transformation, we can then use common "hermitian" perturbation theory to calculate the 
corrections to the energy. 

In fact, in Section|3]a similarity transformation was introduced that does just this, see Eq. 
(I30j) . Although we will take the form of this transformation as a motivation, the transforma- 
tion (J3U|) as it stands is plagued by infrared divergencies. For this reason we rather choose 
a transformation of the Foldy-Wouthuysen type because in this case the Baker-Campell- 
Hausdorff formula guarantees a convenient commutator structure for the transformed oper- 
ator. Since the similarity transform has to map a nonhermitian operator to a hermitian one 
(at least to a certain perturbative order), it can obviously not be unitary. We hence consider 
transformations of the general form 

S = exp (-T (1) - T (2) -...): Q , (99) 

where the operator T (2) is of higher perturbative order than T (1) , etc. It follows for the 
similarity transformed Hamiltonian 

H' B = S-'HbS = H b + [T (1 \ H b ] + \ [T (1 \ [T^,H B ] } + [T (2 \ H B ] + .... (100) 

We could now take the similarity transformation to coincide with the one in Eq. ()30|) to 
lowest nontrivial order in g, 

S> = (U B U B )- 1/2 + 0(g*), (101) 

by choosing 

T' w = \ I dhf dhe-^+^PoH^H^Po (102) 

" J — OO J — oo 

[cf. Eq. (JH2J]- However, this is still a quite complicated choice because of the infrared 
divergencies mentioned before, so we will take a somewhat simpler alternative and leave 
the diagonal contributions (vacuum and kinetic energies) out of T /(1) , defining T (1) explicitly 
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through its matrix elements in the center-of-mass system as 



( P ,r, S |T<vyy> 

[ua(p, r) u A (p', r')} [u B (-p, s) u B (-p', s')] 

2y/2E£2EB2E$2EB E$ + E* - E$ - E* 

x == — = r + 



, , (103) 

E$ + y /(p-p')* + vP-E£ E* + y /(p-p')* + f-Ejj ' 

Definition (|103|) properly applies to Yukawa theory, for the Wick-Cutkosky model one just 
has to omit the products of Dirac spinors. We have introduced a mass \i for the exchanged 
boson as an IR regulator (see below). 

Equation (jlUUj) above, to lowest nontrivial order, then leads to 

H' B = H B + [T« H ] + 0{g A ) = H w + 0(g 4 ) , (104) 

just as for the full transformation in Eq. (J33j) . However, while we were content in Sectional 
with calculating Hy/ strictly to order Hf, for the lowest-order fine and hyperfine structure 
we need all terms of order m r a 4 , and such terms may arise from the order-g 4 contibutions 
to H' B . In fact, we will now show that this is the case, which means that H' B ^ Hw to the 
order m r a 4 . 

In order to extract the terms of order m r a 4 from the order-g 4 terms in Eq. (|100|) . it is 
sufficient to work with T (1) to lowest order in a momentum expansion, explicitly 

(p,r,s\T^\p',r',s') = 2 ™ 5 rr ,5 ss , + 0(a*) (105) 

Kp - p'r + in 

(for Yukawa theory, for the Wick-Cutkosky model just omit the Kronecker deltas). For the 
discussion of IR convergence or divergence, as well as for the estimation of the order in a by 
IR power counting, we have to take into account that Eq. (|105J) is integrated over d 3 p' when 
it is applied to a function of p'. This adds three powers of a to the power counting, so that 
the leading term in Eq. ()105|) is IR divergent and of order a (the IR regulator is ignored in 
the power counting). The fact that T (1) is IR divergent to lowest order is uncomfortable. 
Incidentally, given the relation of T (1) to (U^Ub)' 1 ^ 2 and hence to Uw, this IR divergence 
sheds doubts on the proper existence of the Okubo map. We have not completely analyzed 
the question of the IR convergence or divergence of Uw yet. In any event, the transformed 
Hamiltonian will turn out to be IR convergent to the order m r a 4 considered at present. 
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Now consider the momentum expansion of 
(p,r,s\H B + ^[T w ,H ]\p',r',s') 

I Ijt I / / A 47r " i i 1 P 2 ~P' 2 I r c , , I, 

= (p, r, s\H \p , r , s ) - — ? 1 - == d rr /5 ss / + C(m r Q! 



(p - p') 2 + /i 2 I 8m r v/(p - p') 2 + /x 2 



(106) 



(again omitting the Kronecker deltas for the case of the Wick-Cutkosky model). Observe 
that all the terms are IR convergent by power counting, hence the IR cutoff \i is not stictly 
necessary here. When the commutator with T (1) is taken in Eq. (jlOOj) . the term with Hq on 
the right-hand-side of Eq. (jlOfij) contributes to H\y in Eq. (|1()4|1 . the following term of order 
m r a 2 does not contribute at all up to orders m r a A , and it is the last term of order m r a 3 
which gives the important order-m r o; 4 contribution with matrix elements 

2^ 2 f d*q (p-gHq-pQ 

m r J ( 27 r)3 [(p _ q)2 + ^3/2 [(q _ p , )2 + ^ ^ • (1U7) 

The latter expression is IR finite by power counting, hence the limit \x — * can safely be 
taken. It is then simple to perform the o? 3 g-integration, with the final result 

a 2 1 

{p,r,s\H' B \p',r',s') = (p, r, s\H w \p', r', s') + — 1 6 rr ,5 ss/ + 0(m r a 5 ) (108) 

m r |p — p | 

(without the Kronecker deltas for the Wick-Cutkosky model, as usual). This result coincides 
with the calculation in second-order "antihermitian" perturbation theory, see Eq. (jH^j) . The 
possible physical significance of the "new" m r a 4 -term is discussed at the end of Subsection 
14.31 It appears plausible from the result (J108)) and the consistency of the general formalism 
that this additional term is also generated by the similarity transformation (U b Ub)~ 1 ^ 2 ■ 

To close this appendix, let us have a look at the contributions to order m r a 5 in H' B . It 
is out of the question to present any details of the corresponding lengthy calculations, hence 
we will be content here with quoting the main results. First of all, it can be shown that 
all contributions of this order come from the next-higher terms in the expansions (|105|) and 
(I106|) . inserted in the commutator 



(109) 



The complete result is cumbersome, but we will concentrate here on the IR divergent contri- 
butions that potentially invalidate the similarity transformation. They are explicitly given 
by the antihermitian expression 



47r 2 ct 2 f d?q p 2 — p 12 

(p, r , r', s>)~^r} (^ [(p _ q)2 + , 2] a/ 2[(q _ p02 + /i2] + 



(110) 
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where "~" refers to the equality of the IR divergent terms. It is now possible to remove 
these terms from H' B by introducing a further operator T (2) in the similarity transformation 
HJ), defined by 



57r 2 a 2 f d 3 q 



(p,r,s\T (2) \p',r',s') = / — ^ 5 rr >5 ss/ . (Ill) 

S ' ' ' 3m r J (27r)3[( p _ q )2 + /i 2 ] 3/2 [(q _ p/)2 + //2] y > 

It is conceivable that all possible IR divergencies can be removed in this way, at least re- 
cursively order by order in a. Note, however, that T (1) as well as T (2) are IR divergent 
themselves. 
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